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1 INTRODUCTION 


Retaining wall is a structural element which is intended to support back fill soils, 
e.g. support deep excavations, as abutments of bridge and as anchored retaining 
wall, etc. Retaining walls are known to be simplest structures but there design should 
satisfy some stability requirements. Those are internal and external stability 
conditions. First one talks about structural stability of different parts of retaining wall 
where as second one gives the wall-soil interaction after construction will be in 
equilibrium. While designing a structure, engineers need to adopt many technological 
and managerial decisions at different stages. At present the optimizing design in 
structural elements mainly follows rules which are formed according to experience of 
designer. Most of the designs will be done by assuming cross sectional dimensions 
and material grades and then lead to safety design by compiling them with design 
code requirements. 


Inthis type of designs optimization techniques are best alternative to designs 
which depends on experience. Optimum design is the process of structural synthesis 
which give designs not only safe but economical. In any design optimization is done 
on reliability base or cost minimization. The constraints are formed by using stability 
conditions. Mathematically it is defined as minimization task of constrained function. 
And itis solved by a method based on simple mathematical Program. In recent years 
new developments in numerical methods has given some optimization techniques 
which are conceptually different from traditional methods. These methods are 
termed as non-traditional methods of optimization. These methods work depending 
on biological behavior, molecular swarm of insects and system of neutral biologics. 


Genetic algorithm method is used as a tool in artificial intelligence and computer 
programming. Search is based on Darwin's theory of survival of fittest. It is non- 
traditional search technique in most of cases it is used to find the most economical 
global optimum solutions. This technique starts with a set of design points with 
available variables; it works on the principal of natural genetic and natural selection. 
From the early set of design points, new set of design points are generated and 
weaker points are removed. 
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The genetic algorithm differs from other methods in several aspects: 


% Genetic algorithm works with a coding of the set of feasible point rather than 
the set generated itself. 

% It finds from a random set of points rather than the set itself. 

% Derivatives of the objective function are not used. 

% Ituses random operation in each iteration process. 


2 Genetic Algorithm Method 


2.1 Fundamentals of Genetic Algorithms 


GAs are based on the principle of genetic evolution in which only the best or 
"fittest" individuals of a population survive to reproduce and create the next 
generation. During the process, some genes of the new-born individuals are 
randomly mutated. This operation introduces new changes that could make an 
individual stronger to survive in its environment, and therefore, more likely to 
reproduce and pass those advantages to the next generation. 


In the mathematical analogy a population consist of several design vectors 
(individuals), each design vector is composed of a set of design variables 
(chromosomes) that can be modified within a specific domain or search space 


during the algorithm execution. In the case of single-objective optimization 
problems, a single objective function that assigns a fitness value to the design 
vectors is defined. The natural processes of selection, reproduction and 
mutation are simulated through numerical operations. At every iteration (or 
generation), the best design vectors of the population based on their fitness 
value are selected, combined, and mutated to create a new stronger population. 
The algorithm usually stops after a certain number of generations is completed, 
if no improvement is observed after a specific number of iterations or if another 
convergence criterion is satisfied. 
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2.2 Definitions 


» Fitness Function 


The fitness function determines how fit an individual is (the ability of an individual to 
compete with other individuals). It gives a fitness score to each individual. The 
probability that an individual will be selected for reproduction is basedonits fitness 
score. 


» Constraint Function 


Constraint is a condition of an optimization problem that the solution must satisfy. 
There are several types of constraints —primarily equality constraints, inequality 
constraints, and integer constraints. 


2.3 Genetic Algorithm with Elitism 


The solution to the optimization problem is found by applying a single objective GA 
with elitism. The algorithm used for the present study is briefly described in the 
following steps: 


> Step 1. Initialization 





The algorithm starts at i = O by randomly generating multiple design vectors [e1, e2,. 
.. es] ona search space that is based on their pre-specified domain and ranges. Each 
design vector is denoted as an individual and all together they form the initial 
population Pi. 


» Step 2. Objective and constraint function evaluation 


The objective and the constraint values are computed for all the individuals of the 
current population Pi based on Equations. 


> Step 3. Selection of the fittest. 


The best individuals from the current population are chosen following the 
dominance-based tournament selection operator. This operator takes two random 
individuals from the current population {ea, eb} € Pi to be compared or “fight” 
against each other. 
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The winner of the fight is chosen accordingto the following dominance criteria. If 
both solutions fulfil the constraint, the one with the lowest objective value dominates 
the other. If both solutions violate the equality constraint, the one with the lowest 
constraint value dominates the other. 


After the winner of the first encounter is obtained, a new opponent is chosen 
randomly from the current population to "fight" against the previous winner. This 
battling process is repeated for a specific number of times known as tournament size. 
The winner of the final fight is the champion of the tournament, and the outcome of 
the selection operator. 


After every tournament, the champion is removed from the available pool to avoid 
selecting the same individual multiple times. 


» Step 4. Crossover 


The goal of the crossover operation is to generate new elements by the reproduction 
of the strongest individuals. From the previously created pool Si, two individuals 


Heep ES] are chosen and their decision variables (chromosomes) are combined to 
obtain two new offspring (lex, ey) € Qi. The simulated binary crossover algorithm 
(SBX) is used forthis operation. 





» Step 5. Mutation 


The mutation operation is introduced to the newly created population Qi. The main 
purpose of the mutation is to maintain diversity and avoid the problem of getting 
trapped into a local minimum. The mutation operator produces a small change to the 
decision variables using the polynomial mutation scheme. The mutation is applied 
with a probability of 1/n, so that on average, one variable per individual is mutated. 





Notice that elitism preservation is achieved by passing the selected elements Si 
unchanged to the next population. 


> Step 6. Termination 


The algorithm terminates if the population has converged (does not produce 
offspring which are significantly different from the previous generation). Then it is 
said that the genetic algorithm has provided a set of solutions to our problem. 


Jn 
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Crossover Operator Mutation Operator: 






Select two individuals and swap Select one individual and 







à section of gene between them mutate the genes in it 


Population = Fitness Value ™ Mating Pool = Crossover ® Mutation 


3 Study the guidelines for the design of cantilever 
retaining wall 


3.1 General forces acting on the retaining wall 


The various forces actingon the retaining wall structure are consistent; formulation 
includes active and passive forces onthe front of the toe and base shear key sections 
and the bearing force of the base soil. 


€ Wc= combined weight of all the sections of the R.C.C wall. 


EA 
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% WS zweight of backfill acting on the heel. 
% Wt = weight of soil onthe toe. 
% Q = surcharge load. 


** Pa= force due to the active earth pressure. Pk and Pt =forces due to passive 


earth pressure on the base shear key and front part of the toe section. 
% Pb =force due to the bearing stress of the base shear. 


T2 T4 B9 1] 


Pb 
General forces acting on the retaining wall. 


3.2 Design steps of retaining wall: 


The following analysis and design step has been followed for the writing of MAT LAB 
program. 


LT n 
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> Determination of depth of the foundation: The depth of the foundation 
depends upon the depth at which adequate bearing capacity is available at the 


site. — Do G = ne) 
Lr — y \1l+sing 


Where 


e pois the safe bearing capacity of the soil at the safe at the site. 
e yis the density of the soil at the safe at the site. 
e ( Angle of repose. 


» Height of the wall: The overall height of wall H is equal to height of the backfill 
atthe face of wall plus the depth the foundation. 

> Base width: The width B of the base normally varies from 0.45H to 0.70H, 
depending upon the characteristics of the backfill. It is derived from the 
equation given below, 





2 (Pa Ya t Ps Ys) 
y H 


H = over all height of the retaining wall. 


p,and p, Total horizontal forces on the wall due to active earth pressure. 
y, and y,are the heights of these forces from the base of the footing slab. 


> Stem height: The preliminary stem height is established from overall height of 
the wall H and the base thickness D as the thickness at the base is determined 
from the imposed shear and the moments. 


> Design of the base slab: The heel is subjected to upward soil pressure and 
downward gravity loads due to its self-weight and earth above it. The latter 
being larger causes a hogging moments. The steel placed atthe topface. The 
toe is subjected to the net upward pressure and reinforcement is placed at the 
bottom of the slab. 


^) 
ES 
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» Design of the stem: The stem is designed as a vertical cantilever subjected to a 
triangular or trapezoidal. Main reinforcement is provided on the back fill face of 
the stem. Secondary reinforcement may be curtailed at the top portion of the 
stem. Secondary reinforcement is provided (00.1596 in mild steel and 0.1296 in 
case of HYSD steel. 


3.3 Stability requirement 


> Overturning: To ensure stability against overturning about the edge of the toe 
slab, stabilizing moments due to the gravity loads should be more than the 
overturning moments due to lateral force so asto get a factor so asto get factor 
of safety in the range of 1.5 to 2.0. according to 1S:456 2000, 

» Sliding: The resisting force against sliding should be more than the sliding force 
so as get a factor of safety of 1.55 


EW 
Factor of safety = Po 
Pah 
Where 
YW =Total gravity load. 
Pan = Horizontal components of earth pressure. 
u(=tan@) = friction coef ficient between base and ground. 


@ = angle of repose or angle of internal friction. 


In case the factor of safetyis less than 1.55, a base key may be provided to obtain 
additional resisting force to sliding. The base key develops a passive earth pressure of 


l. 2 
p, = 5 k,yhj;a 


Where h1 is the depth of the base key including the thickness of the toe slab. 


% Base pressure: The vertical pressure on the soil under the base should not exceed 
the bearing capacity of soil. There should be no tension under the base, since 
tension means that. The heel will be lifted above the soil, which is not permissible. 


^) 
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4 Optimization of cantilever retaining wall 


L Optimizer is built on bases of genetic algorithm programming; the problem is 
automated by writing a program in MATLAB software which is user friendly and 
flexible for optimization of cantilever retaining wall. Few important features of the 
program are mentioned bellow. 


% Analysis and designing of the retaining wall is coded in the MATLAB software. 

4 Each commands used in the code are clearly defined in the form of comment for 
self-explanatory of the program. 

*- The user has to enter variables, constraints and the number of individuals per 
subpopulations and maximal Number of generations and generation gap. 

% The user has to enter permissible stresses and permissible deflection for the grade 
of concrete selected as per code. 


L Following are the inputs required for the optimizer: 


% No of variables 

% No of constraints. 

% Maximum number of generations. 
% Lower limit of variable. 

% Upper limit of variables. 


4.1 VARIABLES 


LP 


| 4 x6 
- F i H 
Ti "ua Tui " = 
mE Wr E MEE RT ET 
E ] 
T ai 
s 
XT = x12 -E | | — | 
| XB ! Distance from toe to the thickness of the 
a a a À base shear kev. 
w Tae HEEL x5 ————————————————————————— 
— am Width of the shear key. 
KEY; X11 
x10 
x2 x3 
x9 | | 
x1 Vertical reinforcement in the base shear 


key. 


Variables diagram 
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4.2 Lower and upper limits 
| | Variables  — [| 


X1=Width of the base. 

X2-Toe slab. 

X3-Thickness @ bottom of the stem. 2 (h/0.9)/10 
X4=Thickness @ the top of the stem. 0.2 
X5=Thickness of the base slab. (h/0.9)/10 


X6-WVertical reinforcement in the stem. 


X 7—Horizontal reinforcement in the toe. 0.0012*X5 0.04*X5 
X8=Horizontal reinforcement in the heel. 0.0012*X5 0.04*X5 
"slg ater trom toe to the thickness of the 0.4*h*(12/11) (0.7*h)/0.9 
base shear key. 

X10=Width of the shear key. (h/0.9)/10 
X11=Depth of the shear key. (h/0.9)/10 


X12=Vertical reinforcement in the base shear 0.0012* X10 0.04* X3 


4.3 CONSTRAINTS 


Inequality 
constramts 





— f Mytotal ` 
FOS —| ————_ [= 0 
| Mhtotal | 


rs, -( e forces from passive pressure 
Htotal / 


Moment at bottom of 
Moment at toe 


Ws — Vus <0 
Vt— Vut=0 
(0.12/100)*b*D _At 


(0.12/100)*b*D -As = 0 
(0.12/100)*b*Ds-At = 0 





Ass — (4/100)*b*Ds = 0 





Ast — (4/100)*b*Dt = 0 


Ash — (4/100)*b*Dh = 0 





f (cost) = Cs*Wst + Cc*Vc 
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Nomenclature: 


v/ 


NN S 


V 


Mhtotal = Total horizontal moment of forces that tends to overturn about toe. 
Mvtotal = Total vertical moment of forces that tends to resist overturning about 
toe. 

FSo = Factor of safety against overturning. 

Pmax=Maximum contact pressure at the interface between the wall structure and 
the foundation Soil. Pmin = Minimum contact pressure at the interface between 
the wall structure and the foundation soil. Htotal=Total horizontal driving forces. 
FSs=Factor of safety against sliding. 

B = Base width of the wall. E = Eccentricity of the resultant force. 
Mrs,Mrt,Mrh,Mrk = Flexural strength of stem, toe, heel, shear key. 
Ms,Mt,Mh,Mk= Maximum bending moment at the face of the wall, toe of slab, 
heel of slab, shear key. Vs, Vt, Vh, Vk = design shear carrying capacity at stem, toe, 
heel, shear key. 

Vus, Vut, Vuh, Vuk = Shear capacity of concrete at stem, toe, heel, shear key. 

As, At, Ah, Ak = Area of reinforcement at stem, toe, heel, shear key. ds,dt,dh,dk = 
Effective depth at stem, toe, heel, shear key. 

Cs is the unit cost of steel, Cc is the unit cost of concrete, Wst is the weight of steel 
per unit length of the wall, and Vc is the volume of concrete per unit. 


The below table shows the rates for different grades of concrete and steel, as per 
Hubli- Dharwad schedule rates book. 





gi 


[3 
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4.4 Program for Algorithm 

The program for algorithm to solve the created fitness function involves linking of 
fitness function, defining number of variables has to be defined, stating lower and 
upper bound values, linking of Constraint function and final code for algorithm. The 
program code for algorithm is given below. 





| myfiness.m myConstraints.m | main.m | + 


7 This file can be published to a formatted document. For more information, see the publishing video or help, 
L && Fitness function 


| function y = myFitness(x) 


* y = [x(1) *xi2) ) *[0. 5* (x (3) *x[4)) *($-x(2) ) 7 
4 - end 
5 


Lommand Window 


3,2230 0.5530 0.3000 0.5830 1,2500 


Fitness function 


a) 


jM 
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Bz ley 


D TAB Pl koi Tn AT the e f - al Frictye s wal i 


myfitness m myConstraints.m mainm + 


@) This file can be published to a formatted document. For more information, see the publishing video or help. 





hor 


w c Y (1 e & 
| 


ad 
Ka 
I 


ub 
I 
MonmÉwUuwmpmuornmnouotunm 


n 
l 
n 
D 








èk Constraint function 
function [c,c eq] = myConstraints(x) 


(5-x (2) *x|[2) ) *&0*0.26; 
(20*0.26* |5-x (2) *x(2)0)"2/72); 
3.85*29*x|5) ^2/2; 

(arb-c|*1.5; 

(5-x|2|] -x|2) 42.5) *3g947t |1/6); 

20* (h-x(31) * [52x (2) 2x (2) 42.5) | 
25*x|3| *(5-x121)? 

(25-20) *0.5* [x [4] -x|3) | * ($-x (21); 
Z25*Xx|1]*x|2h:; 


= 20*(xí5)-x(2))*ix|l1l-nl:; 


e* (h-x|[3) } *2.5-£* [h- (x |3) *0. 51 } -3* (h-x [3] - ( [x (4) -X (3) } /3] ) *3*x (1] *0.5-3* [x [32) -h*D. 5) ; 
ky (24-204: -2); 
(e-f-34143)/x(1)7 

k*a* [5-x(2)*x(2) )*D.545* (5-x|2] =x [21] ) *0. 33-c*x|[5) *0.33 se-«f-Cc197|-(x(l1)/2 H 


Command Window 


3.2232 2.5530 0.3000 0.5830 1.250€ 
fval = 
3.7457 


Constramt function 


; >P Users FP sasidnar tanneeru P DESKTOP F Matlab Tes * INESIS P angie OTINternmaltriction P 30 P 3m wall 


A Editor - C:\Users\sasidhar tanneeru\Desktop\matlab files\thesis\angle of internal friction\36\5m wall\main.m 


| myFitness.m = | myConstraints.m =| main.m > | T | 


i» This file can be published to a formatted document. For more information, see the publishing video or help. 





** main code for minimizing fitness function using GA (Genetic Algorithm) 
[= Obj Fen = GmyFirtness: 


| =- nvars = 5; 


| = LB = 
i- UB = 


[3.2 0.5 0.3 0.58 1.25]? 
[3.9 0.562 0.35 0.61 1.269]; 


So ConsFcn = @myConstraints; 
Gis [x,fval] = ga(O0bjFcn,nvarzs,[];,[1;[1,[1],LB,UD,ConsFcn) 


- 


-ommand Window 


3.2230 0.5530 0.3000 0.5830 1.2500 


fval = 


3.7457 


Program for genetic algorithm 


o 


[15 
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9 Numerical example 


A cantilever retaining wall to retain earth for a height of 


horizontal. The density of soil is 18kN/m?. S.B.C of soil is 


of friction between concrete and soil as 0.6. The angle of repose 


concrete and Fe415 steel. 


. The backfill is 
. Theco-efficient 


. Use M20 






1) Parametric studies on effect of change in cost of retaining wall for various grades 
of concrete for 4 mts height. 


Optimum desien variables 


X! | X8 
8 | *410° 


X) | X5 X9 | X10 | XII 


Cost of retaining wall for various grades of concrete 








e Weight of the retaining wall does not show significant changes for 


different grade of concrete. 


e Note: 
% X1, X2, X3, X4, X5, X9, X10, X11, are geometric variable which are 
expressed in meters. 


% X6, X7, X8, X12, are reinforcement variable which are expressed in sq. meters. 


2) Parametric studies on effect of change in retaining wall, heights for M20 grade of 
concrete. 


[16 
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Optimum design variables 


)s 
12 





Optimum design variables for different heights for M20 grade concrete 





e Weight to height ratio is 12.03, 13.013, 18.9, and 18.8. 
e Note: 
% Parameters such as SBC- 200 kN/m2, y= 18 kN/m3, u-0.6, $-30*. 


3) Optimum cost and variables of retaining wall. 


Optimum design variables 


17,350.0 





After optimization of retaining wall 


Optimum cost and variables for 4 m height and M20 grade of concrete. 





e Note: 
** Parameters such asSBC- 200 kN/m2, y= 18 kN/m3, u-0.6, $-30*. 


[Uu 
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Comment : 


** Genetic algorithm can be satisfactorily applied for the optimization of 
cantilever retaining wall. 

* There is sudden change is observed for weight and height ratio, between 
5mts and 6mts 

% Weight of the retaining wall has been reduced from 72.25 KN to 48.13KN. 

* Some geometric parameters of retaining wall have been reduced and some 
are increased. 

* Cost of retaining wall has reduced from 17,350 to 16,000 Rs/mts of length. 


6 CONCLUSION 


Use of genetic algorithm technique in the design of retaining wall has given 
considerable percentage of cost optimization in all the cases which are taken for the 
analysis. This genetic algorithm technique is found to be useful in the designs where 
number of trails have to be done for making it more economical. It requires lower 
and upper bound values of considered variables. The range of boundary conditions 
will alter the objective function value. Hence, the provision of genetic algorithm 
technique in the design of retaining wall is advantageous in optimization of cost. But 
it is time taking method as it requires new program for every change properties of 
retaining wall. There is more scope to work on this modern optimization technique in 
the design of structural elements. 
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